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Abstract 

We calculate the contribution of pions to the (^-expectation value n(p) =< M\qq\M > 
in symmetric nuclear matter. We employ exact pion propagator renormalized by nucleon- 
hole and isobar-hole excitations. Conventional straightforward calculation leads to the 
”pion condensation” at unrealistically small values of densities, causing even earlier restora¬ 
tion of chiral symmetry. This requires a self-consistent approach, consisting in using 
the models, which include direct dependence of in-medium mass values on n(p), e.g. 
the Nambu-Jona-Lasinio-model. We show, that in the self-consistent approach the p- 
dependence of the condensate is described by a smooth curve. The ”pion condensate ” 
point is removed to much higher values of density. The chiral restoration does not take 
place at least while p < 2.8po with po being the saturation value. Validity of our approach 
is limited by possible accumulation of heavier baryons (delta isobars) in the ground state 
of nuclear matter. For the value of effective nucleon mass at the saturation density we 
found m*(po) = 0.6m, consistent with nowadays results of other authors. 


PACS. 13.75.Gx Pion-baryon interactions; 24.80.+y Nuclear tests of fundamental 
interactions and symmetries; 24.85.+p Quarks, gluons, and QCD in nuclei and nuclear 
processes 


1 Introduction 


In this paper we present the calculation of density dependence of the scalar quark condensate 
k =< M\qq\M > in symmetric nuclear matter. The gas approximation [l]] provides linear 
dependence, while account of interactions in medium leads to nonlinear contributions. Although 
at normal density po the nonlinear corrections are rather small Jl|] - 
follow their behaviour while the density increases. 


11 , it is interesting to 
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We show, that the main contribution to nonlinear terms comes from the interaction of 
the scalar quark operator with the pion cloud. Thus we must average the operator qq over 
in-medium pions. 


Concrete calculation are carried out by use of Feynman diagram technique (Fig. 1) with the 
pions described by propagators, renormalized by interactions with medium. The propagation 
of pions in medium is a special story, related closely to the problem of pion condensation |12| , 
T3j. The key point of the latter phenomena is that a solution of dispersion equation with 
negative frequency squared, or? < 0 emerges at certain value p c . This signals the change in the 
structure of the ground state. The general belief is that the phenomena does not take place up 
to the values of about twice the normal density. 


We show that the pion condensation is related to the appearance of a new branch of the 
solutions of pion dispersion equation on the physical sheet of pion frequencies, besides the well- 
known branches with pion quantum numbers (these are the pion, isobar and sound branches). 

Account of the condensation singularity in the calculation of k =< M\qq\M > leads to 
important physical consequence: the value of k turns to zero at the density values smaller than 
p c . Thus, while the density increases, the chiral phase transition takes place earlier, than the 
pion condensation does. 

We show that the value of the critical density p c (and thus, the value of the density of the 
chiral phase transition) depends strongly on the magnitude of the effective nucleon mass m* in 
medium. This happens because the pion polarization operator II is proportional to m*, thus 
entering the dispersion equation and propagator. The effective mass m* decreases while the 
density increases. The decrease of m*(p) causes the diminution of the polarization operator 
reflecting the weakening of the medium influence on the pions. This results in a shift of the 
point of the pion condensation to a larger density. 


However, assuming any of conventional nuclear physics equations for in-medium mass m* 
with direct dependence on the density, we find the chiral restoration point to be dangerously 
close to the saturation value. This would require strong precursors of chiral restoration at 
normal densities, in sharp contradiction to our knowledge. Thus, the intermediate result of our 
paper is that straightforward application of the pion nuclear physics to the calculation of the 
scalar qq expectation value in the nuclear matter leads to unphysical results. The problem is 
solved by self-consistent treatment of the condensate and of hadron parameters. 


Indeed, the expectation value of the quark condensate in nuclear matter is calculated by 
Feynman diagrams technique. The expressions, corresponding to the diagrams, include depen¬ 
dence on a number of in-medium hadron parameters. These are nucleon and pion masses, the 
coupling constant, the number of quark-antiquark pairs in pion, etc. On the other hand, these 
in-medium characteristics can be expressed through the averaged values of quark (and gluon) 
operators in the framework of QCD sum rules []IJ, or by using other models, e.g. Nambu-Jona- 
Lasinio-model (NJL). Such models, combined with the idea of scaling, developed by G.E.Brown 
and M.Rho |14|, [|l^] enable us to express the baryon effective masses m*, m* A and pion decay 
constant /*, through in-medium value of the quark condensate (m* (k), /*(«)). 

Thus, we solve the following system of equations: 


2 





m* = 


with the dots standing for other hadron parameters, depending on k. As a result, there appears 
a self-consistent scheme for the calculation of the expectation value of quark condensate k and 
effective baryon mass m* in the medium. The self-consistent calculation leads to a rapid 
decrease of the effective nucleon mass with density. Thus, there is no pion condensation at 
least up to the density p = 2.8 p 0 . 

At larger densities p > 2.8po the heavier baryons (isobars) can be accumulated in nuclear 
matter, thus changing the structure of its ground state. Therefore, the first phase transition, 
which takes place while the density increases is the condensation of heavier baryons in the 
ground state of nuclear matter. If we neglect the new Fermi sea of isobars, we find that k 
approaches zero asymptotically, i.e. at p —> oo. Also, there is no pion condensation. 

In the self-consistent approach the shape of p-dependence of both functions m*{p ) and n(p) 
does not change much, while we modify the shape of the dependence m*(/t). Dependence of the 
values of m* and k on the values of parameters describing effective particle-hole interactions 
and on those of form factors of the pion vertices is also weak. 


Scalar quark condensate 


Investigation of the scalar condensate can be interesting from several points of view. It may 
appear to be useful in the attempts to find the bridge between description of strong interactions 
in hadron and quark-gluon degrees of freedom. Indeed, the condensate is determined through 
quark degrees of freedom, depending, however on the values of hadron parameters. On the 
other hand, it describes the properties of the matter as a whole. Being the order parameter of 
the system, it characterizes the violation of chiral symmetry. Its turning to zero leads to bright 
consequences for the system as a whole []. 


Also, we hope that the developing of QCD sum rules project started in [||] will provide 
the bridge between the two ways of description, based either on hadronic or quark degrees of 
freedom. By using the QCD sum rules, the particle properties (mass, pole residue, etc.) can be 
expressed through expectation values of the quark and gluon operators. The scalar condensate 
qq is one of the most important. The expectation values of the vector gq 0 g and scalar qq opera¬ 
tors play the role of the vector and scalar boson fields, in terms of the Quantum Hadrodynamics 
(QHD) approach [|T6| in the mean held approximation. In this approximation the nuclear mat¬ 
ter can be considered as a medium with nonzero expectation value, < M\ gq 0 g \M >^ 0, and 
with a new value of scalar operator < M\qq\M >, which is not equal to the vacuum one. 


the other hand, turning of n(p) =< \qq\ > to zero may be not sufficient for the chiral symmetry 
restoration. It is not excluded that the expectation value< M\qq\M > is equal to zero, but simultaneously 

< M | qqqqqq\M >/ 0, etc. However in many of the models used nowadays, e.g. NJL-model, the value < \qq\ > 
may be considered as an order parameter, and the chiral symmetry violation does not take place, although 

< \qq\ > turns to zero. 
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The value of the vector condensate < M\q^q\M >= 3 < M\N ^ a N\M >= 3 p does not 
depend on the choice of degrees of freedom. The latter can be quark or hadronic as well. This 
is due to the vector current conservation. Thus in the latest equality p should be treated as 
the baryon charge density. The density dependence of the scalar quark condensate n{p) = 

< M\qq\M > is more complicated. It is the subject of the present paper. 

Now we run through the main ideas and results of this paper. In the gas approximation [Jp 


k(p) = /c(0) + p < N\qq\N >, 


( 1 ) 


where «(0) =< 0|gg|0 >~ —0.03 GeV 3 is vacuum expectation, and the number of quarks in a 
nucleon, 

, . 2(7 

< N\qq\N >= ---, 

m u + m d 

is given by the well known nN cr-term: <7=45 MeV, m u and m d are the current quark masses. 

— ITJ attempts were made to go beyond the gas approximation 


In a number of papers 


K{p) = K(Q)+p<N\qq\N>+S{p), (2) 

with the nonlinear term S(p). This contribution comes mainly from the pion exchange because 
of small pion mass and large expectation value 

9 777 2 

rj =< 7r\qq\7i >= --—. 

m u + m d 


Being a Goldstone meson, the pion describes the collective mode, i.e. the excitation of the large 
number of gg-pairs. 


The calculations of S(p) in one- and two-loop approximations (i.e. with one- and two-pion 
exchanges) were done in [|I|, |D|, TT] Q At low density the loop expansion is equivalent to the 
expansion with respect to the Fermi momentum pp oc p 1 / 3 . Each extra loop gives an extra 
factor p 1 / 3 (in the case m n -C Pf )■ 


Phase transition 


At large densities, the different types of the phase transitions can take place in nuclear 
matter before the quark-gluon plasma formation. There is a number of possibilities : 

i) new types of baryons (A, £, A) can appear in the ground state of nuclear matter |17|], when 
the energy of the nucleon on the Fermi surface (£f = p 2 p/2m) becomes larger than the mass 
splitting Am = mg — m (B=A, E, A; m = mjv). More precise condition which takes the 
baryon-hole interaction and the mass renormalization into account will be discussed in Sect. 5 
and in a separate paper [fT5|] ; 

ii) chiral invariance will be restored, when a:(p) turns to zero; 

iii) the pion condensation [12] can take place. 

2 Note that erroneous isotopic coefficient was used in Q. So, the result for (pi published in [|ll] should be 

multiplied by 3/4. 
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Let us clarify the latest point. In nuclear medium a pion can be absorbed, producing a free 
nucleon, or A, and a hole. In next step the free baryon can emit a new pion and go back filling 
up the hole. These transitions can be interpreted as the pion interactions with the particle-hole 
channels. Thus, instead of the free pion in vacuum one deals with the mixture of the pion and 
baryon-hole states. These pion- to baryon-hole transitions can be described completely by the 
polarization operator II(u;, k ; p) in the pion propagator D\ 


Here and below k 


D 


1 

uj 2 — m 2 — k 2 — U{co, k ; p ) + is 


\k\. The dispersion equation 


D 1 = J 2 — m 2 — k 2 — n(u;, k; p) — 0 


(4) 


has several solutions. While the density p increases, one of the solutions, co c (p), turns to zero. 
For the first time it happens at a critical density p = p c for some concrete value k = k c of 
pion momentum. At larger densities p > p c , the square of u c is negative ( co 2 (p) < 0 ) in some 
interval of values of k. In this case one has to add the pion-type excitations into the ground 
state of the nuclear matter. This phenomenon is called the pion condensation [12 . 


Since there is a large number of pions in the ground state, large contribution to the expec¬ 
tation value n(p) appears. Hence, the value of k (p) changes near the point p = p c significantly. 
The pion condensation is the main source of nonlinearity in the /t(p) behaviour. 


< \qq\ > in the presence of ”pion condensation” 

The aim of this paper is to calculate the quark condensate in the nuclear medium, with 
account of possible pion condensation. 


We will consider the simplest (one-loop) approximation (Fig. 1) but with the exact (renor¬ 
malized) pion propagator (see Eq. (4)) including geometrical series of baryon-hole insertions. 

Of course, this is not the whole set of Feynman graphs, but it describes and includes all the 
main physical effects we would like to discuss. The short-range interactions will be taken into 
account in terms of the Theory of the Finite Fermi System (TFFS) [|Hj, by using the effective 
constants g' NN , g ' NA , g AA corresponding to nucleon and A-isobar rescatterings. On the other 
hand, the long-range correlations are described by the exact pion propagator. 


Note that in the limit of a small pion-nucleon coupling, our approach reproduces exactly 
H and the most important part of the 2-loop calculations pT . 


the one-loop result 


If the nucleons are treated in nonrelativistic limit, the value of polarization operator, ac¬ 
counting the baryon-hole loop is proportional to the coupling constant g* A / J ’* squared and to 
the nucleon effective mass m*: 


U(uj,k;p) oc m*p F k 2 , (5) 

where m*, g* A and /* are the nucleon mass, axial current and pion coupling in medium, cor¬ 
respondingly. Here and below all effective variables, renormalized in the nuclear medium, are 
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supplied with an asterisk. The Fermi momentum pp oc p 1//3 , and for symmetric nuclear matter 


cpf 


P = 


(2vr) 5 


d 3 k. 


Factor ppm* comes from the integration of the energy denominator, with A E r^j k 2 /2m*: 

rpp d 3 k 


A E 


m pp. 


Thus, to calculate the true value of quark condensate n(p) one has to know the p-dependence 
of the baryon mass m*(p) and that of the coupling constant g\/f*. The simplest possibilities 
are to use either the Landau formula [|20| 


mn 

771 


1 + 


2 Tripp 


fi 


or the Walecka-type model [ 16], where in nonrelativistic limit 


m = m — cp 


( 6 ) 


(7) 


with certain constant coefficients /i and c. 

Using Eqs. (6) or (7), one obtains the pion propagator pole at a > — u ; c (p) with u) 2 < 0 for 
the densities p > p c . The value of critical density p c turns to be of the order of the saturation 
one. 


Say, the value, obtained in the papers |T3|, |2T|, is p c ~ (1.0 — 1.5)p 0 - The appearance of 


such a pole was interpreted as the signal of ”pion condensation’' 0- 

However, just before the ’’condensation” (at p < p c ) the nonlinear contribution S(p) in¬ 
creases drastically, and the curve k( p) crosses the zero line. The reason is trivial. When 
p —> p c , the integral for the pion loop (Fig. 1) takes the form 


S 


dud 3 k 


(u 2 -u ;2(p, k)Y 


(here we keep the singular part of integrand only and take into account the symmetry of pion 
propagator in respect to the sign of u ;). The integral diverges near the condensation point 
a; c (p c , k) ~ a ■ (k — k c ) 2 —> 0 for k —> k c ^ 0: 


duk 2 d{k — k c ) 
(a; 2 — a 2 (k — k c ) 4 ) 2 


oo. 


This means that one faces another phase transition. Namely, the chiral symmetry restoration 
is reached before the pion condensation. At larger densities the pion does not exist any more as 
a collective Goldstone degree of freedom, the baryon mass vanishes (if very small contribution 
of the current quark masses is neglected), and we have to stop our calculations based on the 
selected set of Feynman diagrams (Fig. 1) with exact pion propagator. 
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The structure of pion propagator singularities, various branches of solutions of dispersion 
equation and calculation of the quark condensate n(p) under the assumptions described by Eqs. 
(6) or (7) with the coupling g\jf* = g A / f n = const are described in Sect. 2, 3, 4. 

With conventional values of TFFS constants ( f n , /i, g' NN ,...) we obtain n(p) = 0 at rather 
small densities p ~ (1.1 -j- 1.2)p 0 . This does not look to be realistic. 

Self-consistent approach 

Therefore, we consider another approach, which is a self-consistent one. We obtain the 
expression for quark condensate k which depends on nuclear density, on effective mass m*, on 
the effective constant /*, etc. On the other hand, the effective (renormalized) values of m*, 
f* ... depend on k. 

Recall, for example, that in the framework of QCD sum rules the baryon mass is determined 
mainly by the gg-expectation value. The relation between the mass and n{p) is even more 
straightforward in the NJL-model. The nucleon (and constituent quark) mass is proportional 
to k, and in medium one finds 

m*(p) — G ■ k(p), (8) 

where G is the constant for the four-fermion interaction. Thus, in order to perform self- 
consistent calculations, we have to solve the set of equations 


« = f *(p, 

m*(p) 


IM 

u 


«(0) + p < N\qq\N > +S(p), 


F m {n, p). 


F n (K,p). 


(9) 


In Sect. 5, we use Eqs. (9), together with a hypothesis about g* A / f* behaviour, and calculate 
the expectation n(p). Two types of behaviour of the ratio g* A /f* are considered: 1). g A /fn 


—g A /fn —const , which is the 
f:/U = m*/m |T4 . 


latest version of Brown-Rho scaling [|l^]; and 2). g* A — g A \ 


Here is our main result. For g A /f* —const and 


™*(p) «(p) 

m m ’ m k(0)’ 

we get rather smooth n(p) dependence. The system tries to prevent the chiral symmetry 
restoration at low densities: the ratio n(p)/n( 0) > 0.2 up to p ~ 2.5po ( K (Po)/«(0) = 0.55). 

To study the stability of the results we consider several other possibilities of m*{p ) depen¬ 
dence. These are: 

™*(P) = ( «(p)'\ 1/3 
m \«(0 )) 
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( 10 ) 


and QCD sum rules motivated formula 


m*(p) tz(p) 2.4p 

m k(0) k(0) 


with the second term caused by vector condensate. All the versions with g* A /f*=const show 
the same qualitative behaviour. On the contrary, for g* A = const []and /*/ f n = m*/m we reach 
the critical point n(p) = 0 at a very small value of Fermi momentum Pf — 200 MeV (p < Po/2), 
since in this case the polarization operator increases while |/c| decreases (II oc m */ f* 2 oc 1/k). 
The vanishing of the ^-expectation value at p ~ p 0 /2 does not take place in the Nature. Hence, 
we have to reject this possibility. 

Note that in the present paper we ignore the strange sector, but even in non-strange sector 
new baryons (A-isobars) appear in the ground state of nuclear matter at larger densities. Thus, 
we cannot continue calculations at p > (2.5 — 3.0)po, before the reconstruction of the ground 
state is carried out. 

Since for moderate densities p ~ 2p 0 the effective mass of nucleon rn* becomes comparable 
with the Fermi momentum pp, we account for relativistic kinematics of nucleons. In this 
case we obtained reasonable values of effective nucleon mass m*(po) = 0.6 m and of the scalar 
condensate «(po) = 0.6 k( 0) for normal nuclear density po- 


2 The main equations 

The lowest order contribution to the quark condensate beyond the gas approximation is pro¬ 
vided by the nucleon self-energy graph shown in Fig. 1, with both nucleon and A-isobar 
contributing to the intermediate state. In vacuum all intermediate nucleon momenta P 2 > 0 
are available, but in medium the momenta P2> Pf are allowed only, because of Pauli principle 
(see Fig. 1). Besides, the pion in-medium propagator should be renormalized; this is shown in 
Figs. lb,d by fat wavy line. Therefore, we calculate the contribution of the diagrams, shown 
in Figs. lb,d subtracting analogous vacuum contributions with bare (vacuum) pion propagator 
(Figs. lc,e). 

From formal point of view, n(p) can be calculated as the derivative of the energy density S 
with respect to the current quark mass 

n(p) = d£/dm q . (11) 

The pion-induced part comes from the differentiation of the nucleon self-energy E^r, which 
corresponds to the diagram of Fig. la 

SAP) 

3 As was mentioned in jlj], the low density 
explanation. 


d£, N dm\ 

P dm £ dm q ’ 

evolution from <?a( 0) = 1.25 to <m(po) = 1 has a special 









where the derivative 


d 1 _ 1 

dm 2 (a; 2 — m 2 — A: 2 — Id) (a; 2 — m 2 — k 2 — II) 2 


squares pion propagator, and 


dm 2 

dm n 


mz 


m n 


= V/ 2- 


u q ,,u q 

Here r) is the number of gg-pairs in pion (Eq. (3)). We use linear PCAC equation for the pion 
mass squared m 2 , obtained by Gell-Mann, Oakes and Renner (GMOR) |f23|: 


m 


2 
7r 


< | gg | > ( m u + m d ) 

2/ 2 


( 12 ) 


In the Feynman graph of Fig. 1, the qq operator is shown by the fat black point which i) stands 
for the pion propagator squared, and ii) multiplies the result by the factor rj. 


Note that we average the operator qq over the pion states but not over intermediate baryon 
states. Due to the Ward identity, which corresponds here to the baryon number conservation, all 
the contributions containing < IV|gg|IV > for the nN intermediate state are already accounted 
for by the second term in Eq. (2). 


From technical point of view, this is supported by the following argument. For the nucleon 
in the matter the condensate < N\qq\N > m can be presented as 


< IV|gg|iV > m =< iV m |gg|IV m > + < N\qq\N > 

Here \N > is free nucleon state and \N m > is in-medium nucleon one with energy E] is the 
self-energy. Using the multiplicative character of renormalization 

\N m >=Z l ' 2 \N>, ^ = Z~ 1, 

we find < N\qq\N > m =< iV|gg|IV >, thus proving our assumption. 

For A-baryons (Fig. 1), we deal with the contributions proportional to the difference 

5 =< A|gg|A > — < IV|gg|iV > . 

Basing on the Additive Quark Model (AQM), we assume that 5 = 0. Thus we take into account 
only the pion contribution to fc(p). 

Certainly, in the strong interactions, there is no reasonable parameter for perturbative series. 
We consider the Erst (one-loop) self-energy diagram of Fig. 1 with full (exact) pion propagator 
instead. In other words, we sum up the selected set of Feynman graphs which are responsible 
for the lowest singularity in the p-dependence of n{p). The expression for Sn(p), illustrated by 
Figs. lb,lc is: 


( dE N \ 

V dE ■ 


Sn(p ) 


f d 3 p dud 3 k 
- vSp J (2tt) 3 {2n)H 


(13) 
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(rl NN D 2 (u,k)e{p F - 


^ o(\p-k\-p F ) 

e(p) — u — e(p — k) + id 


-KNND^,k)e(p F -p) 


e(p) — iv — e(p — k ) + id, 


Here Sp stands for summation over the nucleon spin indices, the factor 3 comes from the 
summation over isotopic coefficients and the last term is needed to avoid the double counting 
and subtract the contribution which is already included into < N\qq\N > in the second term 
of Eq. (2), which is related to the bare nucleon in a vacuum. Recall that p = 2p|./37r 2 . We 
must add similar contribution of intermediate isobar, illustrated by the diagram of Fig. Id. 


The free pion propagator is 


D 0 = (a; 2 — k 2 — m 2 + id) 1 . 

The 7T NB vertex with B labeling nucleon or A-isobar is 

RrACB = • ds{k) ■ X-itNBi (14) 


while 


~p(0) 

1 nNN — 




Wa 


X*(vk)x, 


(15) 


with ip(x) being the (non)relativistic nucleon four- (two-) spinors. In order to take into account 
the nonzero baryon sizes, the bare vertex R 0 ' is multiplied by the form factor taken in a simple 
pole form 


1 ~ ml/A% 
' 1 + k 2 /A 2 B 


Aa^O.667 GeV, while Aa=1-0 GeV |24|]. The factors x^nn an d .xy/vA are accounting for the 


renormalization of the corresponding vertices due to the particle-hole pairs. Explicit expressions 
for them will be given later on (see Eq. (27)). 


The expression, corresponding to the diagram of Fig. Id with A-baryon is analogous to Eq. 
(13). However, the unrenormalized vertex is now 


~p(0) 

^ nN A 


Um^=j-x"(Sik)x c 


(16) 


The experiments provide the value of the coupling constant /a/jv — 2 |24|], while AQM calcu¬ 
lations give /a/jv — 1.7. 

Certainly, there are no limits for momenta of isobars in the intermediate states. Mass difference 
Am = ?71 a — m is included into the of A-baryon energy e{p— k). 


We assume that the baryon-medium interactions change the potential energy of any baryon 
(nucleon or A-isobar) by the same value, which does not depend on the baryon momentum p. 
This is consistent with QHD picture in the mean field approximation, under the assumption 
that the vector field has the same coupling to nucleon and A-baryon. In other words, AQM is 
assumed to describe the vector field interaction with baryon. 
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Note that only nonrelativistic approximation for the baryon propagator is used in 
Sects. 2-5. In Sect. 6, to estimate relativistic effects at large densities (when the Fermi 
momentum pp becomes comparable with effective mass m*) we use traditional perturbative 
approach, with all particles on the mass shell. Still, we neglect the contribution from badly 
time-ordered graphs, where the antibaryon-baryon pairs are created. 

Therefore, the only relativistic effect (Sect. 6) is a relativistic expression for e{p) in the 
denominator of Eq. (13): 


£ p = \J (m* 2 +p 2 ). 

Now we construct the pion propagator in nuclear matter. According to TFFS, we have to 
sum up the geometric series of baryon-hole loops shown in Fig. 2a. The contribution of one 
loop, illustrated by first diagrams in Figs. 2b,c are 


n !v = Spj Arjr I kn G n (p + k)9 (|p + k\ - p F )9(p F - p), (17) 

n a =Spj + k)0(p F - p). (18) 

Here the traces are taken over spin and isospin variables, and Gn, Ga are nucleon and A-isobar 
propagators. 


In the described perturbative series for the baryon-hole loop, we consider both particle-hole 
excitation and absorption contributions (the first and the second diagrams, correspondingly, in 
Figs. 2b,c) 


Thus, 


n -- 4 ( 4) 42 


rr(°) _ 
il A — 


<f> 7 v(cc;, k ) + -k) J d? N (k), 

Y ^A/v^ 2 [ $ a(w, k) + $a(-w, -k)] d 2 A (k ), 


with Migdal’s function 


(19) 

( 20 ) 


$a(cu,/c) 


1 m* 3 
47T 2 k 3 



ln( 


a + b. 
a — b 


— ab 


( 21 ) 


Here a = u — k 2 /2m* — Am, b = kpp/m *, and Re(Am) = rri A — m, Im(Am) = —Ta/2. Ta is 
the isobar width. 


Integration over the momenta p provides 


$ N (u,k) 


m* 1 f-um*-\-kpF (kpp) 2 — (com* — k 2 /2) 2 um* — kpF — k 2 /2 
k 47r 2 y 2 2 k 2 um* — kpp + k 2 / 2 


( 22 ) 
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—a nn* ln(- 

um* 


urn* 

— kpF + k 2 /2 


? 


at 0 < k < 2 pp, while 


$jv(w, k) 


m* 1 
k An 2 



at 2 pp < k < oo. 


P/2) + 


(W 2 


(turn* 

2fc 2 


fc 2 / 2) 2 urn* — kp F — k 2 / 2 \ 

um* + kpF — k 2 /2 J 
(23) 


Note that Eq. (22) differs from analogous equation of Refs. 0, 0 at k < 2 pp- The 


short-range correlations originated from the baryon-hole rescattering are described in terms 
of TFFS with the help of effective constants g' NN , g ' NA , g AA , which correspond to N-N, N-A 
and A-A rescatterings. If it is not specially mentioned, we use g' NN =1.0, g' NA =0.2 , g' AA = 0.8 
After summation of the geometrical series of baryon-hole loops, we obtain the polarization 
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operator [T3[], [^3 


n(w, k ; p) = ITat + n A , 


with 


TT° 

11 A 


n JV = n^(i + (7A-7AA) 7 f )/e, 

n A = n° A (i + ( 7A -7 NN )^f)/E. 


Denominator E has the form 


TT° 

11 N 


TT° 

n A 


' JVJV JF2 p2 1 V OVi 

The effective constants 7 are related to g\ as follows: 


)/E, 

(24) 

)/E. 

(25) 

_ 2 \ n ^v n A 

7i) k‘ ■ 

(26) 


7 nn — C 0 g' NN 


'V2 1 ; 


,7a = 


^'Ofl'jVA ( n/2 f* 


_ C 0 g' AA (V2f: 

? TAA p2 


9 a J U/n V 9 a ) fl/N V 9 a J 

where Cq is the normalization factor for the effective particle-hole interaction in the nuclear 
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n- 


C 0 = 

ppm* 

The vertex renormalization factors introduced in Eq. (14) are 

x-kNN — (1 T (7a XnNA — (1 + (7a — 1nn)-^-)/E. (27) 
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3 Singularities of the pion propagator 


We start with the cuts corresponding to singularities of the polarization operator II. In the 
complex cn-plane, the nucleon-hole state induces logarithmic cuts on the real axis. As one can 
see from Eqs. (22-23), U° N has two cuts at positive u in the interval 0 < k < 2 pp 


kp F 


m 


2 ) 

m 2m 

while for large k > 2pp we have one cut only: 

kpF ,k 2 

-b -— < u < 

m 2 m 

The cut caused by the A-hole state lays below the real axis (for Re to > 0) at 


k 2 

2m 

(28) 

kp F k 2 

m 2m ’ 

(29) 

kpF k 2 

m 2m 

(30) 


k 2 

--b Am 

2 m 


kpF , ^k 2 kp F 

-< oj <-b Am H-. 

m 2m m 


(31) 


with Im u> = —T A /2. 

Besides, there are symmetric cuts at Re u < 0. The complete structure of the cuts of 
polarization operator is shown in Fig. 3. 

To make the picture more visual we start with fixed T A =115 MeV equal to the width of 
A-isobar in vacuum. On the other hand, the values of u which are important in onr integrals 
are rather small, and practically there is no phase space for the decay A —> nN in the medium. 
Thus, the A-isobar width in the medium T A ~ 0. Therefore, in final computations we put 
Ia = 0. 

Another set of singularities is provided by the poles of the total pion propagator, i.e. by 
solutions of the dispersion equation 


D 1 = u 2 — k 2 — m* 2 — n(cu, k ; p) = 0. 


(32) 


rn* denotes in-medium value of the pion mass, which is equal to 


m 


*2 

7T 


ml + n, s 


with the scalar polarization operator II S describing S-wave pion-nucleon rescattering. The 
polarization operator, provided by Eqs. (24), (25), accounts for the P-waves only. For the 
analysis, carried ont in this Section it is sufficient to use the gas approximation equation 

< N\qq\N > (m u + m d ) 

n, = -p- M -. 

while in later analysis we use GMOR expression (12) for m* and modify the value of f n as well. 

Starting analysis of Eq. (32) from small values of densities (p < po) we find three branches 
of its solution on the physical sheet. 
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3.1 The pion branch u n (k) 


At k —> 0 it starts at uj — ±yrn* 2 leaving the physical sheet through the cut of 11^. For 
example, this takes place at k ~ 4777^, when p^=290 MeV. Hence, at large momenta one deals 
with the A-hole excitations instead of pion’s ones. 

3.2 The sound branch co s (k) 

It is a slightly changed solution of the equation E = 0 (see Eq. (26)). At very small momenta 
k the admixture of other branches is tiny and the change is negligible. As it should be for the 
sound wave, uj s = const ■ k for small k. At k ~ 0.43777,,. (at p F —290 MeV) this branch leaves 
for the lower sheet through the second cut of U° N (see Eq. (29)). At k > 0.43777,,, this solution 
is on the second unphysical sheet of the complex cn-plane. 

3.3 The isobar branch a ;a (k) 

It is mainly the A-hole sound wave. Starting at Re[uj&(k = 0)) = m A — ra, Im(u)/±(k = 0)) = 
—r A /2, it plunges under the isobar n A cut (at k = 3.8m,,, when p F — 290 MeV). 

3.4 The ” condensate” branch c u c (k) 

This solution comes to the physical sheet through the first cut of at p F > 283 MeV/c. The 
’’trajectories” of the solution are shown in Fig. 4 for different values of p F and of r A . The part 
which is in the upper half-plane (Im uj > 0) corresponds to the unphysical sheet. For small 
Pf A 283 MeV all the trajectory is placed on the unphysical sheet, but at larger p F it comes 
down to the physical sheet. Say, for p^=300 MeV (360 MeV) the solution is on the physical 
sheet at k/m n = 1.06 A- 2.60 ( k/m „ = 0.36 A- 3.91). This is illustrated by Fig. 4b. In Fig. 4a 
one can see that the real part of ui c {k) decreases with T A tending to zero, when u c (k) is on the 
physical sheet. For T A = 0 the solution goes along the negative imaginary axis. Of course, here 
uJc(k) < 0. Thus this is the singularity responsible for the so-called ”pion condensation” 

We have several reasons to put the latest words in the quotation marks. To start with, 
this is not the pion branch (uv), but another one. It starts at k = 0 at the same point as the 
pion branch cov does : u> n (k = 0) = m*. However, it goes to the other sheet. Note also, that 
while Re uj c > 0, the imaginary part of the solution is negative everywhere on the physical 
sheet. This means that we never face the mode with Im uJi > 0. In other words, there is no 
” accumulation of pions”, contrary to the naive understanding of condensation. 

Of course, there are singularities with Im uj t > 0 in the left half-plane of uj. However, these 
singularities originate from the inverse time-ordered graphs of Fig. 2c. They are caused by the 
terms <f>(— u, — k) in Eqs. (19), (20). These singularities correspond to ’’antiparticles” and do 
not describe solutions which grow with time. 

Nevertheless, the fact that even for T A = 0 one obtains a nonzero imaginary part Im uj c ^ 0 
and uo“l < 0 signals on certain instability of the solution. When ujf turns to zero (for any k = k c ) 
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at certain p = p c , the ground state should be reconstructed. New components, like baryon-hole 
excitations (with the pion quantum numbers) emerge in the ground state of nuclear matter. 
Thus, we cannot use the same approach at larger values of p > p c . 

Thus, the appearance of the singularity uR c = 0 on the physical sheet shows, that phase 
transition takes place in the nuclear matter. 

We have to emphasize that in all calculations here and below we use the Landau effective 
mass (6) with the coefficient /i, which gives m*(p 0 ) = 0.8m; then g* A = 1.0 and f* = f n = 
92 MeV. The mass splitting Am = m A —m* = const with Re A = 292 MeV. The values of TFFS 
effective constants are g' NN = 1.0, g' NA = 0.2, g' AA = 0.8. The constant in the 7r!VA-vertex (16) 
is f A/N = 2. 

The dependence of the concrete values of Ui(k) on the values of TFFS constants g' NN , 
g’ NA , and g' AA , being changed in reasonable limits, is weak. We can say the same about the 
dependence on the value of the coupling constant Ja/n- The whole picture is more sensitive 
to the values of effective mass m* and to that pion-baryon coupling g* A jf*. This will be the 
subject of the analysis carried out in next Sections. 


4 Pion contribution to the quark condensate 

In order to carry out integration over c o in the integral in right hand side of Eq. (13), we specify 
the integration contour in complex plane. The contour should go below the pion propagator 
singularities in the left half-plane tv and above the singularities in the right one. We have chosen 
a straight line Im u = a - Re u with a slope a < 1. Since the Cauchy integral is convergent, the 
result of integration does not depend on the slope value a, that is proved by our computations. 

The results of calculation of the function S'(p), defined by Eqs. (2), (13) are shown in 
Fig. 5. For the sake of convenience, we display the ratio 

<P) , . _< N\qq\N > , S(j>) 

MO) I p MO) | MO) I' l ' ) 

The most interesting events take place at pp between 270 MeV and 320 MeV. The large change 
of the values of S(p) is due to the ’’pion condensation” singularity ca c , coming to the physical 
sheet very close to the integration contour at pp — 283 MeV. 

Fig. 5a illustrates the weak dependence of the behaviour of the function S(p) on the values 
of TFFS parameters and on that of Ja/n- 

The value of the width of A-isobar is much more important. At smaller Ta the resonance¬ 
like structure becomes more pronounced (see Fig. 5b). The peak becomes higher and more 
narrow. At Ta —> 0 the poles at u = u> c [k) pinch the contour at u c (k c ) = 0, k c ^ 0 leading 
to the infinite value of S(p), when the ’’pion condensate” singularity emerges for the first time 
on the physical sheet at u c — 0. Recall that namely zero value of Ta is expected in nuclear 
medium for a small uj. 

However, as is was discussed in the Introduction, before the pole at u> = u> c (at p = p c ) 
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reaches the physical sheet, S(p ) becomes so large that it cancels the negative vacuum expec¬ 
tation value k( 0) ~ —0.03 GeV 3 , and the whole scalar quark condensate turns to zero. The 
vanishing of the scalar quark condensate indicates the chiral invariance restoration, and after 
that one has to deal with quite another system, i.e. with another phase of the nuclear matter. 

The prediction of chiral phase transition at rather low values of p (close to normal nuclear 
density p 0 ) looks too strong. On the other hand, the results are stable enough and do not 
change too much under the variation of TFFS parameters and that of /a/tv- These statements 
are true for both types of behaviour of nucleon effective masses. However, the dependence on 
the numerical values of the coefficient which enter Eqs. (6), (7) is strong. The coefficient can be 
fixed by the choice of the value of m*(po). The results for m*(po) = 0.8m and for m*(po) = 0.7m 
are compared in Fig. 5c. In the latter case, for the smaller effective mass, the phase transition 
takes place at pp — 320 MeV, i.e. at larger values of density. 


5 Self-consistent approach 

5.1 Assumptions on the density behaviour of hadron parameters 

As we showed above, the phase transition (either ”pion condensation” or chiral symmetry 
restoration) density value depends strongly on the value of baryon effective mass m*. On the 
other hand, in the framework of commonly used models the hadron mass depends mainly on 
the scalar quark condensate ft(p). This problem should be solved self-consistently. 

Of course, it would be nice to calculate all the masses (m*, m* A , m*) and constants ( g\, 
/*,..) with the help of QCD sum rules, substituting them in the next step into our expression 
for k(p), solving the equation 


k = /<p,m*(p, A),/;(p, A),...) (33) 

in the final step. Unfortunately, it is not so easy. One of the main obstacles is that the masses 
and hadron constants depend not on the value of the scalar quark condensate only but on the 
in-medium expectation values of other operators (usually, more complicated ones) as well. 

Therefore, we use simplified scenario of the k dependence of the mass m* and of the other 
parameters involved. 

Fortunately, there are model-independent equations. The GMOR relation can be generalized 
to the case of finite density due to PCAC: 

m ** = ~ ( m « + m d)■ (34) 

Also, due to chirality, the number of quark-antiquark pairs inside pion is 


2 mf 

(m u + m d ) ‘ 


(35) 
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Now we come to model-dependent relations. Denote the ratios 


In the NJL-modcl 


11 


a(p) = 


m 


m 


P(P) = 


k 

f: 


( 36 ) 


a(p) 


«(P) 

«(«)' 


(37) 


In order to check the stability of our result, we perform the calculations using the two more 
types of ^-dependence of a : 


a(p) 


kp)\ 1/3 

<0)J 


(38) 


This latter expression can be justihed by the dimensional counting, if there is only one dimen¬ 
sional parameter. Unfortunately, in our case, this argument does not work, since there are at 
least two external dimensional parameters p and A qcd, 


Another expression for a(p) is motivated by 




a(p) 


«(p) _ 2A P 
k(0) k(0) 


(39) 


with the last term in the right hand side caused by the vector condensate. 

Now we discuss the in-medium behaviour of pion decay constant /*, which is proportional 
to the pion radius inverted, i.e. f n oc a jN~ c /r n [^6[. In the NJL-model near the phase transition 
point (k —> 0), the pion radius increases unlimitedly, and f n —> 0 when k =< \qq\ >—> 0. Thus, 
it looks natural to assume that 


/tt m ' 


Brown and Rho JlJ| made even stronger hypothesis assuming that all the parameters of the 
same dimension are proportional to each other 


rrS = m^ = fl 
m m A //" 


a(p), i.e. (3 = -, 
a 


(40) 


while the dimensionless parameters do not change in nuclear medium; in particular, (see foot¬ 
note 3) 

9*a = 9a = const. 


The second, alternative hypothesis is based on the idea of confinement. If the deconfinement 
phase transition does not take place simultaneously with the chiral invariance restoration, the 
pion radius should be limited and f* has a nonzero value when we approach the chiral transition 
point. Therefore, we consider below another limiting possibility: 
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fl = U, i.e. p = 1. 


(41) 


In the present calculations we have fixed the value of the axial coupling constant, g* A = 1. The 
most important parameter in our calculation is the ratio g* A /f* (see Eq. (5), (19), (20)). Thus, 
just as in the latest version of Brown-Rho scaling []nj, g* A /f* = const, i.e. does not depend on 
P- 


Note that both hypotheses, described by Eqs. (40), (41) are consistent with the QCD sum 
rules for pion ]27|], which can be generalized for the case of finite density in a straightforward 
way (the proof will be published elsewhere): 


1 ( f>l V 3wy / q,(wy) \ 8/l ’ n a, 2 

2 \m u +m d ) 32tt l a.Qi) ) 16 ir [ 1 

Here Wq is the continuous threshold value, i.e. the minimal energy of the multihadronic states 
with pion quantum numbers; a s is the QCD coupling and G 2 is the gluon field squared. 

Neglecting the last (numerically small) term and anomalous dimension (i.e. putting C( s (Wq) 
—a s (p)), one can satisfy Eq. (42) in two ways: (i) f* — — const and the threshold position 

Wq 2 oc m* 2 , or (ii) the fixed threshold W 0 * — Wo — const 1 GeV and Eq. (40) is consistent 
with for f* oc n(p) 

We must also make assumptions on in-medium value of the A-isobar mass. If Eq. (40) is 
true, the A-isobar-nucleon mass splitting satisfies the relation 

(m* A - m*) _ Am* _ _ 1 

(tua — m) Am (3 


However, if /3 = 1 (Eq. (41)) we come to Am* = Am. 

The experimental situation with A-isobar mass in nuclear matter is not quite clear at the 
moment. On the one hand, the total photon-nucleus cross section indicates that the mass m* A 
does not decrease in the medium |^S|, while the nucleon mass m*(p) diminishes with p. This 
means that the splitting Am* increases (Am* > Am), opposite to the p — n mass splitting. 
This fact (if it does take place) looks strange, since the two kinds of splitting are caused by the 
same colour magnetic (spin-spin) quark-quark interaction. On the other hand, the experimental 
data for total pion-nucleus cross sections ^9f are consistent with the mass m* A decreasing in 
the matter. As to calculations, the description within the Skirmion model |JD| predicts that 


m* A decreases in nuclear matter and Am* < Am. Equation Am* = Am is also true in Walecka 
model, if AQM prediction for the scalar field-baryon coupling g S NN = g s aa is assumed. 


Thus we expect that the phenomenological parametrization, 


Am* _ 1 

Am /3(pY 


i.e. Am* = Am when (3 — 1 does not look too unlike. 
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5.2 A rejected scenario 


Here we try the scaling, provided by Eq. (40), with a(p) given by (37) and (3 = 1/a. We find 
the phase transition to take place at the densities, about 2.5 times smaller than the normal 
one, (pf — 200 MeV). 

The technical explanation is simple. Polarization operator given by Eqs. (5,19,20,24-26), 
which is responsible for the ”pion condensation” singularity ui c (k), behaves as 


n 


m 

^ 7*2 ^ 

J 7T 


1 

«(P) 


oc 


1 

a (p)' 


While the density p increases, the value of |/c(p)| and that of a(p) become smaller. Polariza¬ 
tion operator increases and the ”pion condensate” singularity c o c approaches the physical sheet 
at smaller p F . The nonlinear pion contribution S(p) becomes very large and the whole value 
k(p) tends to zero (dashed curves in Fig. 7). Also, since n oc l/a(p) —> oo when the value 
of k{p) turns to zero, we lose the solution of equation (33) for n(p) (see Appendix for details), 
which from now on becomes the complex one. 

Since that, one has to deal with another phase of the matter, with much smaller particle 
masses. The whole picture contradicts sharply to our knowledge about nuclei and nuclear 
matter. Thus, we reject this scenario. 


5.3 An accepted scenario 

For (3 — 1 (Eq. (41)) the situation looks much better. Up to p ~ 2.5po, we deal with a self- 
consistent solution of the system (9). The value of quark condensate tends to zero, when that 
of the density p increases, and we never reach either the chiral symmetry restoration or ”pion 
condensation”. 

Self-consistent set of Eqs. (9), is solved numerically using the standard iteration procedure. 
Technically this means, that calculation of n(p) with vacuum parameters is followed by calcula¬ 
tion of a(p). In the next step we obtain in-medium values of all the other parameters by using 
Eq. (40). This enables us to obtain and Ha provided by Eqs. (24)-(27). In the last of the 
cycle substitution of these operators into Eq. (13) for further integration provides a new value 
of k(p). Thus we come to a new cycle. 

The validity of our calculations for the larger values of p is limited by another phase transi¬ 
tion. At some value of Fermi momentum (p FA oc p 1 / 3 ), die total energy of a nucleon on Fermi 
surface becomes larger than the energy of the A-isobar at rest, i.e. E A (0) < E N (p F ). Thus, 
the isobars starts to be accumulated by the ground state of nuclear matter. This effect, as 
well as possible appearance of other types of baryons, can be taken into account in our scheme. 
Thus our calculations are reliable below p F , corresponding to this phase transition and mark 
the corresponding points by the black circles on the curves shown in Figs. 6, 7. 

From the technical point of view, the effect of instability with respect to the A-isobar 
accumulation in the ground state reveals itself in the fact that the branching points (left edge 
of the right A-hole cut (with Im u < 0) and right edge of the left A-holc cut (see Fig. 3)) 
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cross the vertical axis Im cu=0. And two A-hole cuts start to deform, i.e. to cross or to pinch 
(for r A =0) the integration contour in the cu-plane (Eq. (13)). 


The possibility for a other, but nucleons, types of baryons to be contained in the nuclear- 


matter ground state was discussed 
condensate was considered in 
in terms of Walecka model in 


m 


r?H . The case of A-isobar in the presence of a pion 


311, [32]; the problem without the 7r-condensation was studied 


The results plotted in Fig. 6a (in terms of n{p)) and Fig. 6b (in terms of m*(p )) do not 
change too much under reasonable variations of the TFFS couplings. Instead of g' NN = 1-0 in 
the master version (solid curve), in Fig. 6 we put g' NN — 0.7 (dot-dashed curve); for dotted curve 
g AA =1.2 instead of 0.8; for dashed curve /a/iv=1-7 instead of 2.0. In all calculations in this 
Section we have used T A =0. 


The lower solid line in Fig. 6b (and Fig. 7b) is drawn for the limit values of m*/m at every 
value Pf■ If the ratio m* jm is smaller than the limit value, the ground state of nuclear matter 
contains isobars. The equation for this line is determined by the condition that the isobar 
logarithmic cut (Eq. (31)) of the polarization operator I1 A (Eq. (20)) touches the vertical axis 
Im oj = 0 see (Fig. 3). 

In Fig. 7 we demonstrate the dependence of the value of the scalar condensate on the type 
of the scaling function a(p). Solid curve corresponds to the master version, with cc(p) given by 
Eq. (37) and (3 — 1. Two other parametrizations, provided by Eqs. (38), (39) shown by dotted 
and dot-dashed curves. The results, corresponding to the law, described by Eq. (38), differ 
quantitatively from the two others. However, as we said earlier, the latter are better based. 


6 Account of relativistic kinematics 

At moderate densities of about twice the normal value, the value of effective mass m* becomes 
comparable with that of Fermi momentum. Thus, one cannot neglect the relativistic effects 
any more. We take into account relativistic kinematics by using relativistic expression for the 
energies = yj ( m* 2 + k 2 ) in all the formulae. However, we still omit the baryon-antibaryon 
pair contributions. 

In terms of traditional perturbative theory, Eqs. (21)-(23) should be replaced now by: 


®(u>,k) = <E>jv (u j,k)9(p F - k) + ^\ou,k)6(2p F - k)6(k - p F ) + (uj,k)9(k - 2 p F )\ (43) 


m CPf CPf r k ~PF 

d’jy (c o,k) = / dpA(tv,k), (ca, fc) = / dpA(uj, k) + / dpB(co,k), 

JpF—k Jk—pF JO 

^(a;, k) = I dpB(u, k), 

Jo 
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where 




A{uj, k) = 


B(u>, k ) = 


/■? 

= 

Jo 


p m* 

* 

m _ 

47t 2 k 

(p 2 + m* 2 ) 

p m* 

m* 

47t 2 k 

yj (p 2 + m* 2 ) 

p m* 

(m* + m* A ) , 


)y (p|, + m* 2 ) — u; — »/ (p 2 + m* 2 ) 


'((p + A:) 2 + m* 2 ) — a; — y (p 2 + m* 2 ) 

y ((p — fc) 2 + m* 2 ) — a; — ^(p 2 + m* 2 )\ 
j{{jp + k) 2 + m* 2 ) — lu — J (p 2 + m* 2 ) / 

((p — k ) 2 + m A ) — uj — y f (p 2 + m* 2 ) 
4 712 ^ 2 \J(p 2 + m* 2 ) \\J{{p + k) 2 + m^) — a > — \J{p 2 + m* 2 ) 


(44) 


(45) 


(46) 


By using these relativistic expression, we can extend the self-consistent calculation up to p ~ 
2.8po- The limit is still determined by transition to the isobar accumulation phase. The results 
of calculation of the scalar condensate for the three considered possibilities of dependence of the 
effective mass on k are presented in Fig. 8. In the QCD sum rules motivated parametrization, 
described by Eq. (39) with (3 = 1, we find m*/m = 0.67, at normal nuclear density p = po- 
And we have m*/m = 0.6 in parametrization (37) with (3=1. The value is in good agreement 
with the one, obtained recently in framework of QHD |35[. 


7 Summary 

Since the gas approximation equation for the quark condensate n(p) was presented in |1]|, the 
nonlinear contribution S(p) was considered in a number of papers. The analysis of M.Ericson 
et al. 0. i was based on the general properties of the ttN scattering amplitude and its 
generalization for the case of nuclear medium. As usually, some uncertainties come from the 
fact that one deals with the off-mass-shell amplitude. Thus, certain assumptions about the 
AW-interaction and on p-dependence of the effective pion mass m* in medium are needed. 

Another group of papers |J, 0, 0 was based on NJL-modcl. However, in this case they 
discussed not the nuclear (build up of the hadrons) medium but the quark one (the quark 
plasma). In such approach the pion constant /* tends to zero, while the mass m* —> oo near 
the point of chiral invariance restoration. As we discuss in Sect. 5, this looks unlikely for the 
real nuclear matter. 

Our approach is based on using the exact pion propagator, renormalized in nuclear medium 
by the insertions of the nucleon-hole and A-hole loops. The short-range correlations are ac¬ 
counted for by methods of TFFS. The lowest laying singularity in p corresponding to the 
so-called ’’pion condensation” is included. We carried out self-consistent calculations with the 
quark condensate n(p) depending on the effective mass m*(p), while the mass m*(p) itself is 
determined by (or strongly depends on) /?(p). Nonlinear p-dependence of k is obtained by 
calculation of the diagram, shown in the graph of Fig. 1. Calculations include the in-medium 
values of nucleon, isobar and pions masses and other parameters (/*, g% g'..). On the other 
hand, QCD sum rules and NJL-model give the relations between masses and quark condensates 
which can be used to determine the in-medium masses and parameters, if n(p) is known. This 
enabled us to solve the set of equations (9). 
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It should be emphasized that, since the dependence of m*(p , n(p)) on n(p) was treated self- 
consistently, the ”pion condensation” singularity was pushed out from the physical sheet. The 
only effect which can limit the validity of our calculations at large densities is the accumulation 
of isobars in the ground state of nuclear matter at p ~ 2.8 p 0 . In the general case one should 
include the possible accumulation of hyperons. To understand, which of the condensates ap¬ 
pears earlier, one should have better knowledge of hyperon interactions with matter. Analysis 
of the problem was started by Pandharipande 0 - This goes beyond the scope of our paper. 

At small densities the nonlinear term S(p) diminishes the absolute value of k(p) in compar¬ 
ison with the gas approximation ( the tendency which was already noted in [|1TJ for very low 
p). In our self-consistent approach such behaviour continues up to p ~ 2.0po- 

Taking relativistic kinematics into account, we have calculated the expectation values of 
scalar quark condensate in the symmetric nuclear matter up to p ~ 2.8po; where k reaches 
the value ft(2.8po) — 0.1 k( 0). At normal nuclear density p = po = 0.17 fm -3 , we obtain 
k(Po) — 0.6/c(0) and the effective nucleon mass m*(po ) — 0.6m. The latter result is very close 
to the value used nowadays in QHD |35| to describe the properties of nuclei. 
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7.1 Appendix 

To clarify what happens in the case of scaling, described by Eq. (40) we omit dependence of 
polarization operator on pion momentum k using the mean value of polarization operator II 
instead. The quark condensate |«(p)| obtained in such a way is plotted in Fig. 9 as a function 
of II for two values of density. One can see the dependence to be a monotonous one. The curves 
Al and B 1, calculated for p = Pa and for p = p B > Pa are shown as Al and B1 in Fig. 9. 
At some value II = II C , when the ”pion condensation” singularity u c comes the physical sheet, 
the function k(p, II c ) —> oo. While both second and third terms of Eq. (2) increase with p (for 
S(p) an extra multiplicative factor p comes from the integral over the nucleon momenta), the 
function n(p, II) becomes steeper for a larger density. To find self-consistent solution in this 
simplified approximation, we have to solve the set of equations 


k = k(p, n), (47) 

n = n(p,«). (48) 

We show the function 11 ( 10 , k) by curves A 2 and B 2 in the same figure. The crossing points of 
Al and A 2 (of B 1 and B2) provide the solution. 
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It was shown in Sect. 1.2 that the polarization operator can be approximated as n(ft,p) ~ 
p 1//3 / n. So II —> oo at n —> 0 and the mean value of II increases with p. 

At low densities, there are two solutions: point 1 and point 2, where the curves A 1 and A2 
cross. For p —> 0 the first one (Ai in point 1) matches smoothly with the solution of NJL gap 
equation in vacuum: K\ —> «(0) when p —> 0. When the density p increases the solutions 1 and 
2 draw nearer, merge with each other (at some p = p c ) and go out into the complex plane. In 
this way we lose the real solution at p > p c . 

Assuming that the baryon radius (i.e. parameter A in the form factors d n, (see Eq. (14)) 
follows the same scaling law ( A*/ A = f* / f n = ...), one obtains more complicated function II(/c) 
shown in Fig. 9b. Now for a very small k the form factors djsr,dA cut off the integral in the 
right hand side of Eq. (13) at k ~ A << pp. Therefore, the effective polarization operator 
II(/c) reaches its maximum value and falls down when k —> 0, proving us with the third solution 
with a very small quark condensate |/c| = k 3 and thus with very small values of hadron masses. 
This solution does not disappear even at very large densities p. On the other hand, for the 
third solution somewhere at a — k%/k( 0) < 0.25 the pion mass 


becomes larger the nucleon mass 


m 


*2 



«(Q)) 

K 3 


1 

OC — 
«3 


l»(0)| 

and we cannot any more apply our approach, based on the summation of a selected set of 
Feynman diagrams corresponding to the lightest pion degrees of freedom. 

In any case, at p > p c we lose the primary real solution 1 and face a first order phase 
transition (mass of a hadron has a discontinuity) which looks rather strange. 
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Figure captions 


Fig. 1. a) Diagrammatic representation of the interaction of the operator qq (fat point) 
with the pion held. Straight line denotes the nucleon; the wavy line stands for the pion. 
b,c) Diagrammatic expression for Eq. (13) with nucleon in the intermediate state. Thick wavy 
line denotes the pion propagator renormalized due to baryon-hole excitations in the framework 
of TFFS. d,e) Diagrammatic expression for Eq. (13) with A-isobar (double solid line) in the 
intermediate state. 


Fig. 2. a) Full pion propagator in the medium (thick wavy line) equal to the sum of the 
geometrical series of the nucleon-hole and isobar-hole excitations, b) Polarization operator of 
the pion, fl^, consists of two terms corresponding to excitation and absorption of the nucleon- 
hole pair, c) Pion polarization operator, 11^, consists of two terms corresponding to excitation 
and absorption of the isobar-hole pair. Thin wavy line denotes a free pion, solid line with left 
arrow is a hole, solid line with right arrow denotes a nucleon, double line is A-isobar. 


Fig. 3. Positions of singularities of the polarization operator fl in the cu-plane (represented 
at pf = pfo and k = m n ). On the right half-plane tu/m n , in the interval w/m, r = 0.0 4- 0.26, 
there is the Erst logarithmic cut of the fl^-function, see Eq. (28). Moving after thin arrow 
from the physical sheet across the first cut, we continue the movement on the upper logarithmic 
sheet (thick arrow). In the interval u/m n = 0.26 4- 0.45 there is the second logarithmic cut of 
fljy, see Eq. (29). In this case, following the thin arrow across the second cut we come to the 
lower logarithmic sheet (dashed arrow). Logarithmic cut of the n^-function is located at 
Re cjj/m tt = 1.824-2.54 and Im u = —Fa/2 = —0.115/2 GeV. Since U° N and 11^ are symmetrical 
in lv i-y —uj permutation (Eqs. (19), (20)), there are symmetrical cuts on the left half-plane 

Fig. 4. The condensate branch u c (k)/m 7V . a) u c (k) is presented at pp = Pfo for isobar 
widths values Ta= 0.001, 0.01, 0.05, 0.115 GeV (curves 1,2,3,4, correspondingly). Solid lines 
are disposed on the upper (unphysical) sheet of the first logarithmic cut. Continuation of 
this branch to the physical sheet is displayed by dashed curves. All curves begin at k = 0 in the 
point u = m* = O.SOm^ (Eq. (32)). b) c o c (k) is shown at r A = 0.115 GeV for different values 
of Fermi momenta Pf=280, 290, 300, 360 MeV (curves 1, 2, 3, 4, correspondingly). For curves 
3 and 4, the part of uj c {k) on the physical sheet is drawn only. The curve 1 for p j p=280MeV is 
completely on the unphysical sheet. 

Fig. 5. The function S'/|/c(0)|. a) Dependence of S'/1 ac( 0) | on the variation of nuclear 
parameters. Solid curve represents the main result obtained with g' NN = 1.0, g\ VA = 0.2, 
g' AA = 0-8, /ttTva = 2.0, Ta = 0.115 GeV; the other parameters are described at the end of 
Sect. 3. Dashed curve corresponds to the calculation with f^NA = 1-7, dotted curve to g' AA =1.2 } 
dot-dashed curve to g' NN = 0.7. b) Dependence of Sy|/s;(0)| on the isobar width. Solid curve 
represents the main result (Ta = 0.115 GeV). Dotted curve corresponds to the calculation 
with Ta = 0.07 GeV, dot-dashed curve to Ta = 0.05 GeV, dashed curve to Ta = 0.01 GeV. 


25 



c) Dependence of S'/|k( 0)| on the behaviour type of m*(p). Solid curve stands for the result 
with m* provided by Landau equation (6) (■ m*(p = po) = 0.8 m). Dashed curve corresponds to 
Walecka equation (7) for m* (■ m*(p = po) = 0.8 m). Dot-dashed curve is obtained in framework 
of Walecka model, with m*(p = p 0 ) = 0.7 m. 


Fig. 6. Self-consistent results for the qq expectation value /c(p)/|/c(0)| =< M\qq\M > /|/c(0)| 
(Eq. (2.1)) and m*(p)/m in nuclear matter, a) Dependence of /c(p)/|/c(0)| on the variation 
of nuclear parameters. Scaling functions are (3 — l,a — k(p)/k( 0). Solid curve is the main 
result obtained with g' NN = 1.0, g' NA = 0.2, g AA = 0.8, f n NA = 2.0. Dashed curve corresponds 
to the calculation with f^NA—lA^ dotted one to g AA = 1.2, dot-dashed curve to g' NN = 0.7. 
b) Dependence of m*(p)/m on the variation of nuclear parameters. Notation of curves are the 
same as in Fig. 6a. Straight line is drawn for the limit values of m*/m (see the end of Sect. 
5). Corresponding reliability limits for the calculation of k/k{ 0) are marked by black points. 


Fig. 7. Self-consistent results for k ( p )/| k ( 0 )| and m*(p)/m in the nuclear matter, a) De¬ 
pendence of tt(p)/|tt(0) on the type of scaling functions a(p),(3(p). Solid curve displays the 
main result obtained with (3— 1, a is calculated by Eq. (37). Long-dashed curve corresponds 
to the gas approximation given by Eq. (1). For dotted curve (3 = l,a is taken from Eq. 
(38). For dot-dashed curve [3 — l,a is taken from Eq. (39). Dashed curve corresponds to the 
Brown-Rho scaling: (3 — 1/a, a — k(p)/k( 0) (40). b) Dependence of m*(p)/m on the kind of 
scaling functions. Notations of the curves are the same as in Fig. 7a. Straight line is drawn for 
the limit values of m*(p)/m (see the end of Sect. 5). Corresponding reliability limits for the 
calculation of k/k{ 0) are marked by black points. 


Fig. 8. Self-consistent results for /c(p)/|/c(0)| and m*(p)/m with relativistic corrections, 
a) Dashed curve for /c(p)/|/c(0)| is calculated with relativistic corrections (Sect. 6) and (3 — 1, 
and a = k/k{ 0), Eq. (37). Dotted curve: (3 = 1, and a is taken from Eq. (38). Dot- 
dashed curve: (3=1, and a taken from Eq. (39). Solid curve represents the main result 
((3 = 1, a = k/k( 0), Eq. (37)) without relativistic corrections, b) The results for m*(p)/m 
calculated with relativistic corrections. Notations are the same as in Fig. 8a. Straight line 
restricts acceptable values of m*/m from below, it is obtained using the same method as 
described above (see the end of Sect. 5) but with the isobar polarization operator, Eq. (46). 
Corresponding reliability limits for the calculation of k/k{ 0) are marked by black points. 


Fig. 9. The graphical solution of the system of equations (47) and (48). a) The curves Al 
and A2 are plotted for the equations (47) and (48), correspondingly, at certain small density 
Pa- The crossing points 1 and 2 are the two solutions of a system of equations. Dashed curves 
B1 and B2 correspond to the density value ps > Pa- There is no solution in the latter case, 
b) Graphical solution of the set of equations (47) and (48), when baryon radius obeys the scaling 
equation (40). Solid curves Al and A2 are plotted for Eqs. (47) and (48), correspondingly, at 
small pa- There are three crossing points for solid lines, which are the three solutions of the 
system. Curves B 1 and B2 are calculated at p = ps , Pb > Pa- In this case only one solution 
survives. 
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